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Abstract 

We argue three-flavour neutrino mixing. We consider the neutrinos as Majorana 
particles and see how the neutrinoless double beta decay constrains the neutrino 
mixing angles. Our formulation is widely valid and is applied to the neutrino os- 
cillation experiment. 
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It is one of the most important problems in particle physics whether the neutrinos 
have masses or not. From the recent neutrino experiments |l| |Q] Q @j it becomes 
very probable that the neutrinos have masses. However, if the neutrinos have masses, 
we must explain the reason why they are so tiny relative to the charged lepton masses. 
Seesaw mechanism is one of the most promising candidates for such an explanation. In 
this case, neutrinos become Majorana particles. 

In this paper we consider the neutrinos as massive Majorana particles with three 
generations and see how this point of view constrains physics of lepton sector. 

As is well known, the neutrino oscillation does not distinguish Majorana neutrinos 
from Dirac ones. So, lets us first consider the neutrinoless double beta decay ((/3/3)ov) 
which occur only in the case of Majorana neutrinos. 

Fig.l 

The decay rate of {(3(3)oi, is, in the absence of right-handed couplings, proportional to 
the "averaged" mass defined by |5| 

3 

(m^) = l^U^jUijl. (1) 
i=i 

Here Uaj is the left-handed neutrino mixing matrix which combines the weak eigenstate 
(a = e, ^ and r) to the mass eigenstate with mass rrij (j=l,2 and 3). It takes the 
following form in the case of Majorana neutrinos, 

/ C1C3 sicse'f S3e^(^-^) \ 

U = (-S1C2 - ciS2S3e'^)e-'^ cica - siS2S3e''t> 82036'^'"'^'^ . (2) 

V {siS2 - ciC2S3e'^)e-'P (-C1S2 - siC2S3e^'^)e-^(^-'3) C2C3 / 

Here Cj = cosOj, Sj = sinOj {9i = ^12, 62 = O23, 63 = 631), and beside (j) , appear the 
two additional CP violating phases /? and p for Majorana particle. So {m,^) becomes 

(m,) = \micl4 - m2s\cle-^'P' - m3sle-^''>' I (3) 

where we have introduced 

p'^\-{p-<t^). (4) 

CP violation occurs in the presence of imaginary part in (mjy), though this process 
itself does not explicitly show CP violation. 

We show that the neutrino mixing angles are constrained from the presence of CP 
violation. Here we follow the method given in Q. 

From Eq.(3) it follows that 

{rriv)'^ = (micf c| — m2s\c^cos2[)' — m3s\cos2p')^ + {m2s\c^sin2j3' + m3s|sin2p')^ (5) 
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Rewriting cos2p' and sin2p' by tanp' , we can consider Eq.(5) as an equation of tanp' , 

a^ptan^p' + bptanp' + = 0. (6) 
Here and bp are defined by 

— A ■ 2 pi 2 2/ 22i 2\|/ 22 22i 2\2 / \2 

a±j3 = 4stn p m2SiC^[miCiC^ ± mss^) + [rriiCiC^ — rn2SiC^ ± m^s^) — {rrii,) {1) 



So, the discriminant D for Eq.(6) must satisfy the following inequality: 

D = bp- Aa^pa^p 

= 4?{m^c{cl)\m2slclf{U - sir? l3'){siT? (3' - /_) > 0, (8) 



where 



So we obtain 



{{my) ± m^j slf - {micjcl - m2sf ci)^ 
4mim2cfsfc3 



/± = A „2,2^4 ■ W 



/_ < sir?P' < /+. (10) 

It follows from Eq.(lO) that 

/_<!, /+>0. (11) 

Quite analogously, rewriting cos2(3' and sin2(3' by tanP', and considering Eq.(5) as an 
equation of tan^', we obtain other inequalities, 

5_ < sir?p' < g+. (12) 

Here 



9± = A „2„2^2 • \^^) 



_ {{vrii,) zb m2sf Cg)^ - (micf C3 - rn^s'l)'^ 
So we get 

5-<l, 5+>0. (14) 



The conditions (11) and (14) are consistency conditions. CP violating area is given by 
the more stringent condition 

< /- < sir?(3' < /+ < 1 or < g- < sir? (J < 5+ < 1. (15) 

From the inequalities (11) and (14), we obtain the allowed region of the mixing angles 
in the s\ versus S3 plane once the neutrino masses mi and the "averaged " neutrino 
mass {my) are known. The magnitude of {m^) is experimentally unknown at present. 
The neutrino masses may be safely ordered as mi < m2 < mz- So in the following 
discussions we consider the three cases: 

a) {my) < mi, 

b) mi < {my) < m2. 
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and 

c) 7712 < {rriy) < 7713. 

Note that the definition of {m^) in Eq.(l) and the Schwartz inequality leads us to 

3 3 
{rriy) < mj\Ulj\ < ms ^ jC/^l = rns, (16) 

so {nil) can not be larger than m^. The allowed regions in the s\ versus 53 plane for 
each case (a), (b) and (c) are obtained from Eqs.(ll) and (14), and are shown in Fig.2. 

Fig.2 

From Fig.2, we obtain the upper bound on S3 as 

si < a^il^ (17) 

7713 + m2 

for any case. The CP violating areas in the sf versus S3 plane given by Eq.(15) are also 
indicated by the oblique lines in Fig.2 for each case (a), (b) and (c). The above case 
(a) was considered also in and[^. In the representation for the mixing matrix 
adopted by Cabibbo-Kobayashi-Maskawa was used. In Q, only the limiting case where 
all the neutrino masses are degenerate (mi = m2 = m^) was discussed. It should be 
noted that we consider the cases (b) and (c) in addition to (a) and that no condition 
on the neutrino masses has been imposed so far. 

The above mentioned method is not restricted to Majorana particles but is widely 
applicable. Next, we consider the constraint from the neutrino oscillation experiment 
at CHORUS Q and see how this method also gives the allowed region in sf versus S3 
plane. We assume here 6m1-^ = mg — ~ (Jm^x = 771-2 ~ "^1 ^ <^"^32 = '^i ~ ^2- 

In this case the approximate oscillation probability is given byP] 



P(z.^ ^ ur) = 4\U,,\'\Uri\'sinH^-^). (18) 



Substituting the expression of Eq.(2) into Eq.(18), we obtain the following equation 
w.r.t. cosd>. 



A = 



a+ cos^ (/) — 26 cos (/> + a_ (19) 
/(cos (/)). 



Here 



a+ = -(2siS2S3CiC2)^, 

= i44 + 4slsl){slsl + cl4sl), (20) 

b = SiS2S3CiC2(Si - Cis|)(c2 - S2). 
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The oscillation process does not distinguish Majorana neutrino from Dirac one, and 
only (p phase takes place. Firstly the discriminant D of Eq(19) leads to 



Thus we obtain 



D = b'^-a+{a--A) 

= isiS2S3CiC2f[{sl + clslf-iA]>0. (21) 



> A, (22) 



which is irrelevant to 023- Therefore the CHORUS data on P{i^^ Vt) (= PcHORUs) 
constrain the allowed region in the versus S3 plane (Fig. 3). 



Fig.3 

Unfortunately, we have only the upper bound on the PcHORUS-, PcHORUS < 2.5 x 10 
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[10|. Setting L = 600m (the midpoint of the maximum length, 800m and the minimum 
length, 400m), E = 27GeV and 5ml^ ~ dml^ = GeV^ » 6ml2, we have A < 0.022. 
The broken line is the trajectory of 

= 0.022, (23) 

and the allowed region is the upper part from the broken line. If the PcHORUS gives the 
lower value, the broken line moves downward to extend the allowed region. The shaded 
areas in Fig.3 are those of Fig. 2(b) under the assumption that Sm'^i ~ ^^-21 = Se^^ ^ 
(5m32 and mi -C (m,y) with possible {m^) values. The more stringent constraints, 
though they depend on 623, are also obtained from Eq.(19). a+ is negative definite 
and /(±1) are positive definite. Therefore from the condition that —1 < cos(p < 1 we 
obtain the following inequalities: 

Case a-1 : < ^^a-'^DC^f-'^j^i) < i 

— 4S1S2S3C1C2 — 

(siC2 - c,S2Ssf{siS2 + C1C2S3)' < ^^''^ ^^^^ ^ I^^^ + ^^^s)' (24) 
Case a-2 : 1< ^''2-4)i^'-^>l) 

^ 4siS2S3ClC2 

{siC2-CiS2S3f{siS2 + CiC2S3f < ^^^^ < {siC2 + CiS2Ssf {siS2 - 010283^ 

(25) 



Case b-1 : -1 < (^2 '^2)(^i '^i^.O < q 

— 4siS2S3CiC2 — 

(siC2 + ciS2Ssf{siS2 " ciC2Ssf < ^^""^ ^J^^^ ^ 7(^? + ^^^D' (26) 

4sin2(^L) 4 
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Case b-2 



(,2_,2)(,2_,2,2) 



4siS2S3ClC2 



^ < -1 



(S1C2 + C1S2S3)^(S1S2-C1C2S3)^ < 



Pil'u Vr) ^2/ \2 

^^-r^ < (S1C2-C1S2S3) (S1S2 + C1C2S3) 

4sin2(?^L) 



(27) 



As we have mentioned, we have experimentaUy only the upper bound of P{i'^ 



at present. So the more meaningful inequalities than Eq.(22) comes from the lower 
bounds of Eqs.(24) ~ (27). Namely we have 



Case a 



(S1C2 - ClS2Szf{siS2 + < 



Case b 



(S1C2 + CiS2S3)^(siS2 - CiC2S3)^ < 



4sin^(5|^L) 



4sin^(?*L) 



/or isi-4)isi-cisi)>0 (28) 



/or (s^-c^)(st-cfs^) <0 (29) 



From these inequalities (28) and (29), we obtain another allowed region in the sf versus 
S3 plane for a fixed value of 62- Using 5m'^i ~ <Jw-2i = 6eV^ » (^'^32 and A < 0.022, 
we show the allowed regions for ^2 = 0, ^, j^, • • • , f in Fig.4. 

Fig.4 

Lastly we comment on another neutrino less process of fi~-e'^ conversion [11]. (Fig.5) 

Fig.5 

In this case the averaged neutrino mass which will be determined experimentally is 
given by 



{mu)n-e+ = 



fll 



i=l 



i=l 



Substituting the expression of Eq.(2) into Eq.(30), we obtain 

{mu)i^-g+ = ml{c1clclsl + 2 cos 0cfc2C3.s1.s2.s3 + cfc3.S2.S3) 



(30) 



I A( Z A A A n ± A 6 I A 4: A A\ 

+m2(c^C2C3S]^ — 2 cos (pCiC2C^SiS2Ss + C^SiS2S^) 

, 2 2 2 2 
+m3C3S2S3 

+mim2(— 2cos(2/3)cic|c|si — 2cos(2/3 — (/')cf C2c|siS2S3 
+2cos(2/3 + 0)ciC2c|sfs2S3 + 2cos(2/3)cic|sfs|s|) 

+mim3(— 2cos(2p — (^)ciC2C3.si.S2S3 — 2cos(2p — 2(j))c1c^S2sf) 

+m2ms{2cos{2P -2p + 0)ciC2c|siS2S3 

-2cos(2/3-2/9 + 2(^)cis?s^si). (31) 
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In contrast to the neutrinoless double beta decay, all the mixing angles and the phase 
parameters appear. So if we assume one of the phases and 623, we can develop the 
same argument as that in the neutrinoless double beta decay. 

In conclusion, we have proposed the new method to constrain the neutrino mixing 
angles from the observed data of {rriu). Our method, however, is widely valid and 
have been applied to the neutrino oscillation, having given the new constraints from 
the observed data of P{iyfx — > i^r)- Our method will be applied to the other decay and 
oscillation processes. 
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Figure Captions 
Fig.l: Feynman diagram of neutrinoless double beta decay. 

Fig.2: The allowed region on sin^9i2 versus sin^ ^31 plane by the neutrinoless double 
beta decay is given by the shaded areas in the respective case: 

(a) {nil,) < mi (b) mi < {m^) < 777-2 (c) 7772 < {rrii^) < rris 

In the allowed region, CP-violating area is specially indicated by the oblique lines. 

Fig.3: Each allowed region by the neutrinoless double beta decay in the respective 
case: 

(777^) = 0.4, 0.8, 1.2, 1.6, 2.0, 2AeV 

under the assumption that Sm^i ~ (^Tiii = QeV^ » <^'T732 and 7771 -C (777,^). 
Broken line is given by Eq.(23). The upper part from the broken line is allowed 
by the neutrino oscillation experiment at CHORUS. 

Fig.4: The allowed regions by the inequalities of Eqs.(28) and (29) under the condition 
that PcHORUS < 2.5 X 10~^, 5m1i ~ 5m2i = 6eF^ S> 5ml2 with given O2. 

(a)^?2=0,f (b)02 = iL,J^ (c)02=-,5. (d)^2 = f,f 
(6)^2 = 1,1 (f)^2 = lf,|| (g)^2 = f 

Fig. 5: Feynman diagram of the pi~ — conversion. 
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->- - 

(A,Z) (A.Z+2) 

Fig.l 
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1112 + nil + rm 




m2 - mi m2 + rni 




Fig.2 
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Fig.4 
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(A,Z+2) 



(A,Z) 



Fig. 5 
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